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Abstract 

We study the trapped surface produced by an off-center coUision of hght-hke, point- 
sourced shock waves in anti-de Sitter space. We find an analytic expression for the shape 
of the trapped surface in the hmit where the energy of the shock waves is large and the 
impact parameter is not too large. We use the area of the trapped surface to estimate a 
lower bound on the entropy produced in the collision. We compare our results to particle 
multiplicity measurements in heavy-ion collisions as interpreted through the Glauber model. 
In an attempt to roughly simulate the effects of asymptotic freedom and confinement in 
quantum chromodynamics, we also consider the effects of slicing off parts of anti-de Sitter 
space. 
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1 Introduction 

The Relativistic Heavy Ion Collider (RHIC) collides gold nuclei at a/s/vtv = 200 GeV. This 
means that each nucleus has energy E = 100 GeV per nucleon, for a total of about 19.7 TeV 
for each nucleus. The total number of charged particles N^harged that emerge from such a 
collision can be as large as 5000: see for example [1]. In |2] we pointed out that this number 
can be approximately reproduced starting from collisions of gravitational shock waves [S] [U 
[5l[6l[71[8l[9l[T0l[IIl[T2]in AdSs, following the methods of the gauge-string duality [131 [H [I5] 
and identifying the total energy of each nucleus with the energy of the corresponding shock 
wave. The calculations of |2] rely upon finding a marginally trapped surface in the five- 
dimensional geometry and using its area to put a lower bound on the entropy of the black 
hole produced from the collision, following the general plan of |16l [TT] [181 [E]- In [2], we 
translated this entropy bound into an approximate lower bound on A^charged using a fairly 
well established relation]^ 

S ~ T.SA'charged 5 (1) 

^The simplest justification for this relation is that it holds, approximately, for a thermally equilibrated 
gas of non-interacting hadrons at temperatures just below the transition temperature Tc « 170 MeV of QCD 

[2ni[ni[ia. 
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between Ai'charged and the entropy S produced in a heavy-ion colhsion 
A virtue of the trapped surface calculation is that there are no free parameters: the gauge 
coupling doesn't enter to leading order in a strong coupling expansion, and the overall 
normalization of the entropy is fixed by the equation of state for static plasmas (see for 
example [2S])- On the downside, the predicted lower bound on Nchavged scales as s]^f^, which 
is a faster energy dependence than the s]('^ scaling predicted by the Landau model p6] and 
largely obeyed by the data. Since the trapped surface computation gives a lower bound on 
the entropy produced, there is no confiict between p| and experiment — so far. Confiict will 

1 /3 

arise if the growth of Ai'charged remains slower than s significantly above ^/snn = 200 GeV. 

In this paper we have two main aims. First, we want to generalize the methods of [2J to 
handle off-center collisions of gravitational shock waves in AdSs. Our generalization amounts 
to giving an approximate answer to a well-posed question: What is the area of the marginally 
trapped surface lying on the past light-cone of the shock waves? The answer, for trapped 
surfaces that are significantly biggeij^than the radius of AdSs, is 

where z=k is a characteristic transverse length scale, and P = b/2z^: is an AdSs version of 
the impact parameter. A more general version of ^ appears as ( [56| ) (see also (58)) and 
represents our main analytic result. 

With ^ in hand, the standard entropy estimate for the black hole created in the collision 

is 

Using ([T]) and fixing parameters in the same way as in [2j , we can obtain approximate lower 
bounds on A^^chargcd and compare them with heavy-ion data — or, more precisely, to heavy-ion 
data as interpreted using the Glauber model. We find that our bound on A'charged has a signif- 
icantly weaker dependence on the impact parameter than what data plus Glauber indicates 

1 /3 

for quantum chromodynamics (QCD). This result, coupled with the s^^y dependence of the 
bound on Ai'charged from trapped surface calculations, may indicate that the agreement found 
in ^ was to some extent fortuitous — or that the best motivated mapping between quantities 
in QCD and quantities in a strongly coupled field theory is more subtle than we proposed 



^More precisely, a sufficient condition for equation ([2]) to be a good approximation is that C = 
{2Ez,^f'^ > 1 and /3< C- 
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in [2]. In fact, we find a substantially improved fit to the data by identifying the energy 
of each shock wave with the fraction of the energy of the nucleus carried by nucleons that 
participate in the collision]^ Figure [s] shows how our total multiplicity estimates compare 
with data from PHOBOS [17] with the original identification of energy as the total energy 
of the nucleus. Figure |6] shows the improved fit obtained from the rescaled energy. 

Our second aim is to inquire how our calculations might change if we took into account 
the non-conformal nature of QCD. It's probably because QCD confines that we have to 
exclude the energy of the non-participating nucleons before we get good agreement between 
AdS/CFT calculations and total multiplicity data. We are led to ask, is there a less con- 
trived way of including the effects of confinement? Low-energy processes in QCD should 
not contribute as much entropy as in a conformal theory: the reason is that the number of 
degrees of freedom decreases dramatically below the confinement transition. To incorporate 
this in our trapped surface calculation, we discard the part of the trapped surface below a 
certain fixed depth, corresponding to an infrared cutoff. Also, for sufficiently high energies, 
asymptotic freedom dictates that interactions become weaker and weaker. So it seems safe 
to say that there is not as much entropy production from hard processes in QCD as there 
is in a strongly coupled conformal theory. Correspondingly, we slice away the part of the 
trapped surface which is above a fixed depth, corresponding to an ultraviolet cutoff. In 
summary: we go back to identifying the energy of the shock wave with the total energy of 
the nucleus; we find the trapped surface in pure AdSs without any cutoffs; and finally, we 
slice away both the ultraviolet and infrared parts of the trapped surface before converting 
its area to an estimate of the entropy produced in the collision. 

This slicing approach is distinct from a hard-wall construction of a holographic dual of 
QCD in that the shape of the trapped surface doesn't respond in any way to the cutoffs. 
This approach is admittedly naive. Our main defense of it is to note that it is still more naive 
to neglect violations of conformal invariance altogether. In a proper treatment, we should 
replace the infrared cutoff with a holographic renormalization group flow to a confining 
theory whose equation of state matches that of QCD: see for example [271 l2H] for work 
along these lines. The ultraviolet cutoff should in principle be replaced with some hybrid 
description where Einstein gravity rolls over into perturbative quantum field theory as one 
approaches the boundary of AdSs. 

1 /6 

The slicing approach described in the previous paragraph leads to a scaling of the 
lower bound on Ncha.Tged at large energies. Heavy ion collisions at the LHC may be able to 

■^We thank B. Cole, J. Noronha, P. Steinberg, and B. Zajc for suggesting this approach. 
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probe this scaling. 

Our analysis is hardly the only attempt at relating heavy-ion collisions to black hole 
formation from collisions of gravitational shock waves in AdSs. Other work along these lines 
includes [HI |29l |30l [311 |2l |32l [33] and references therein. One hazard of our calculation is that 
black hole formation is complicated, and it's not clear whether the inequality S > S'trappcd is 
close to being saturated. For collisions in flat space, the recent study [3l] suggests that it is 
not: instead one has S > 1.5 S'trapped, where the inequality is approximately saturated in the 
limit of ultra-relativistic collisions, and S'trapped is based on the canonical choice of trapped 
surface, analogous to the one we use. Another hazard is that heavy-ion collisions are a multi- 
scale process involving both perturbative and non-perturbative dynamics. A supergravity 
treatment of an AdSs dual can at best give us a reliable handle only on the range of energy 
scales where QCD is deconfined but still strongly coupled. Slicing away parts of AdSs is a 
poor man's approach to dealing with the non-conformal nature of QCD. 

The rest of this paper is organized as follows. In section [2] we describe the shock wave 
solutions we need and exhibit a helpful 0(2) symmetry of the collision. In section [3] we 
find the shape of the trapped surface in a limit where the size of the surface is much larger 
than the radius of AdSs. In section |4] we use the optical Glauber method to relate our 
gravity computations to the experimental results of [35]. In section [s] we describe in detail 
our approach of slicing off the UV and IR parts of the trapped surface to obtain a revised 
estimate of A^^charged that accounts, albeit crudely, for the non-conformal behavior of QCD. 
We end with a discussion of our results in section [6] 

When this paper was nearing completion, we received [36], which has some overlap with 
our results. We will comment on the relation between our results and the ones in [36] at the 
end of section |3} We also received [37] , which takes a somewhat different approach to shock 
wave collisions in AdSs, following the earlier work [H2] . 



2 Shock waves from point particles 

Point-like massless particles moving in AdSs generate gravitational shock waves: they are 
sources for the geometry in the sense that the five-dimensional stress-energy tensor has delta- 
function support on null geodesies. From the boundary point of view, a collision of two such 
particles will correspond to a collision of localized energetic objects which we take to be 



nuclei. In section 2.1 we discuss the shock wave geometries before the collision and explain 



their field theory interpretation. Most of the material in this section is more thoroughly 
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reviewed in [2], where references to the original hterature are also given. In section 2.2 we 
introduce coordinate systems that make the symmetries of the problem manifest. We show, 
among other things, how conformal symmetry relates off-center collisions to collisions of 
objects of unequal size. 



2.1 Shock wave geometries 

Particles in AdSs can be described by the action 

s ' 





12" 




1 d^x^ 




+ J dr] 



dx'i dx'i e 9 
m 



2e^''" df] dr] 2 



(4) 



where x^ = x'^{r]) is the trajectory of the particle, m is its mass, r] is an arbitrary parametriza- 
tion of the worldline, and e is the worldline one-bein. We of course set m = 0. Ignoring the 
back-reaction of the particle, the geometry is pure AdSs: 



+ {dx^ + {dx'^f + {dx^f + dz^] 



(5) 



and the particle's trajectory is a null geodesic in this geometry. The momentum conjugate 
to xt is 

(6) 



1 dxl 
Pfi = -g^lu- 



e ' dr] 

Because of the translation symmetry in the x'" directions, where m runs from to 3, the 
quantities Pm are conserved. They can be identified as the fiat space four-momentum in the 
gauge theory. Let's focus on a particle whose trajectory is 



x^^ = xt = t xi = x1 = Q z^ = constant 



(7) 



Then we have 

= {-E,0,0,E), 

where E is the energy of the particle. Introducing the light-cone coordinates 

x^ = x^ ± x^ p± = ^{po± ps) , 



(9) 



one finds from ([8j) that p_ = —E. We may choose r] = t as the worldline parameter. Then 
from ( 6 ) one finds that = 
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Now let's add in the back- reaction. The Einstein equations following from Q are 



where 



J^^ = I dr] -^=5^{x^' - x'i{r]))p^p^ 



(10) 



'11^ 



is the stress tensor of the particle. Using light-cone coordinates, the only non-zero component 



of J^,y is 



Plugging the ansatz 



J__ = E—6{x^)6{x^)6{z - z,)6{x') 



ds =dsj^^g A $(x ,x ,z)(5(x ){dx ) 



:i2i 



(13) 



into (10), one finds that the only non-trivial equation is the equation, which reads 

^H, -J:^^ = -167rG5i?|^5(xi)5(x2)5(z - z,) , 



(14) 



where 

v2 = \(—\ (A 

~ L2 \dx^ J ^ \dx^ 
is the laplacian on H^, whose line element is 

L2 



+ 



dz^ 



z d 

T^lfz 



(15) 



ds]j.^ 



[{dx^ + {dx^f + dz^] . 



(16) 



The differential equation (14) is subject to the boundary condition that $ as one 



approaches the boundary of H^. This is equivalent to requiring the perturbation of the 
metric to vanish at the boundary, meaning that the shock wave describes a state in the dual 
field theory rather than a deformation of the lagrangian. 

Reducing the full non-linear Einstein equations in five dimensions to a single linear differ- 
ential equation on ifs is a drastic simplification. It is important to realize that no linearized 



approximation is needed to derive (14): a solution to it leads to an exact solution of the 



Einstein equations. Only, because 6{x ) is involved in the metric (13), it is a solution in the 
sense of distributions. 



The solution to ( 14 ) is 



G^Ez, 1 

8L2 g3 



,Fi(3,5/2;5;-l/g) 



(17) 
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where 



and 2Fi{a, b; c; z) is the hypergeometric function. The particular hypergeometric function in 



(17) has a closed form expression which is algebraic in q, but its explicit form is not very 
illuminating. The geodesic distance (in H3) from the point (0,0, z*) to where the point- 
source is located is i = 2Lsinh~^^/g. Thus, the solution to (14) is invariant under 0(3). 



This 0(3) is a subgroup of the isometry group of H^, which is 0(3, 1)/Z2|^ 

One may extract an expectation value for the gauge theory stress tensor dual to the point 
source by using a standard expression for one-point functions: 

= 7^S{x-) hm ^i^^llp^ = l^'* 5{x-) , (19) 
47rG5 z-^ n[Xj_ + z^y 

where x± = {x^, x^) parameterizes the plane transverse to the collision. All other components 



of {Tmn) vanish. One can straightforwardly check from (19) that 



/ 



£x{T__) = E and ^ f = ■ (20) 



The first equality shows that E is the total energy of one shock wave in the gauge theory, 
as it is in the gravitational description. The second equality shows that the energy-weighted 



root-mean-square size of the distribution (19) is z^:. This is an elementary example of the 



relation between depth in AdSs and size in the dual field theory. 

Now consider an off-center collision of two shocks, where the trajectories of the point 
sources take the form 

= ^x° x^ = b± x^ = z = z±. (21) 
A straightforward choice would be to set b± = ±6/2 and z+ = z_ = 2;=,,: then we would be 



colliding an energy distribution like (19) with another one of the same size going the other 
way. But we can just as easily consider z^ z^, which means that we are colliding energy 
distributions with different transverse sizes. 

Causality dictates that the gravitational shocks generated by these particles cannot affect 
one another outside the future light-cone of the collision point. So one must be able to 



^0(3,1) has four connected components. In the hyperboloid coordinates, to be introduced in (32 1, the 
Z2 we divide by acts as — > —X~^. Dividing by an additional Z2 acting as parity on the remaining 
coordinates would give 50(3, 1). 
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superpose them to obtain an exact (distributional) solution of the Einstein equations which 
holds provided x+ < or a;^ < 0: 



By the same calculations that led to (17), one finds 



$^ = ^^^^4 2Fi(3, 5/2; 5; -l/g±) . (23) 

Here i?+ and are the energies of the two shock waves, which we do not assume to be 
equal. The quantities q± are the chordal distances between a given point {x^,x^,z) and the 
location {b±,0,z±) of the shocks: 

= -l^^ • ^2^) 



Passing through the calculations that led to (19), one finds that the gauge-theory shock wave 
moving in the direction generates a boundary theory stress-energy tensor 



77 [(Xl - 6_)2 + (a;2)2 + ^2 ]^ 



(^--) = . , .i3 ^(^-) > (25) 



while the shock moving in the —x direction generates 

2^4 , , , 

TT [(X^ — 0+)2 + (x2)2 + zW 

While it is difficult to see it in the Poincare coordinate system, the two-shock geometry 



(22) has an 0(2) symmetry. In the following subsection we use a different coordinate system 
which makes this symmetry explicit and will also make it easier to deal with collisions of 
unequal-sized energy distributions in the boundary theory. 



2.2 Symmetries of the collision 

In ([5]) and ([o]), we have introduced the Poincare patch coordinates for AdSs and our conven- 
tions for light-cone coordinates. However, the 0(4,2) conformal symmetry becomes more 
transparent if one uses the hyperboloid coordinates X^^ in R^'^, which are subject to the 
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constraint 



- {X-'Y - {X'Y + {X'Y + (X^)' + {X'Y + (X^)' = -L\ (27) 

By using 

X™ = L— {21 



2 V " ■ z"^ 

we see that the metric inherited on this hypersurface from the standard flat metric on R^'^ 
is the same as ([s])]^ 

The most general collision we will consider is that of two shock waves whose point sources 
follow the trajectories 

X^ = TX^ X^ = ±Lf3cosa X^ = X* = ±L(3sma X'^ = L^l + (3^. 

(29) 

Here /3 is a dimensionless AdSs version of the impact parameter. Changing a evidently 



amounts to a conformal transformation in the gauge theory. In Poincare coordinates, (29) 
becomes 



where 



= ^x" = h± = z = (30) 



LB cos a L 

b± = ± , z± = . . (31) 

^l + /?2±/5sina ^/lT]3^±psma ^ ' 



By choosing a = 0, we find ourselves colliding objects of equal size; but if a 7^ 0, the 
objects have different sizes. In short: the relative sizes of colliding objects can be changed 
by a conformal transformation! Experiment mostly focuses on colliding equal-sized nuclei, 
but it would be interesting to inquire to what extent the collision of unequal-sized objects 
(say, gold against copper) would provide experimental tests of the degree to which conformal 
symmetry is preserved in the dynamics of a heavy-ion collision. The residual 0(2) symmetry 
which we mentioned earlier manifests itself as rotations which leave the quadratic form 
[X^ sin a — X^ cos a)^ + (X^)^ invariant. 

Most of our analysis will take place on i^s, which is the intersection of the null surfaces 



^It's worth noting that global AdSs is usually thought of as the covering space of the hyperboloid (27 1, 
which has closed timelike curves. 
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= and x = in the Poincare patch. From (28) we see that is determined by 
= X^ = — almost. If we impose these two conditions on (27), we obtain the equation 



-{x-y + {xy + {xy + {X 



-4\2 



(32) 



which describes a two-sheeted, three-dimensional hyperboloid in R ■ . But (27) also implies 



that X~^ > in the Poincare patch, and this condition selects the upper sheet. Sometimes H3 
is denoted to emphasize this restriction, but we will instead use to mean just the upper 
sheet. Setting = = in (28) leads immediately to a coordinate transformation between 



hyperboloid coordinates (X ^,X^,X^,X^) on and Poincare coordinates {x^,x'^,z). 
To make the residual 0(2) symmetry of transparent, we switch to radial coordinates 



X^ = r cos 9 cos a — r sin 9 cos sin a 
X"^ = r sin 9 sin 

X"^ = r cos 9 sin a + r sin 9 cos d) cos a . 



(33a) 
(33b) 
(33c) 
(33d) 



in which case the line element on takes the form 



-dr^ + r\d9^ + sin^ 



In this coordinate system the chordal distance (24) takes the form 



cos 9 

= ~2 + 2L ^ ~2L~ 



(34) 



(35) 



and the 0(2) symmetry of the solution is made clear by the invariance of (23) under rotations 



of 0. (Recall that setting a = corresponds in gauge theory to colliding objects of equal 
size, separated in the x^ direction.) 



3 Trapped surface computation 

If the impact parameter is not too large, then after the collision a black hole will probably 
form. To obtain a lower bound on the entropy of the black hole produced in such a collision, 
we use the method of [16] , which was further developed in [TTl [TH [191 EH] and adapted to 
AdS space in [5l EJ [TJ HI [TOl [391 121 [10] . The method is to find a marginally trapped surface 
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iS, composed of two parts: 5 = 5+ U 5_ where the iSj's are parameterized by 



5± = {x'' e AdSs : x"^ = 0,x^ = -^±(a;^x^2)} , (36) 

and then use the area theorem and the Cosmic Censorship conjecture to give a lower bound 
on the entropy produced: 

Q \ C — ^trapped 

•J ^ ^trapped — ^ ^ 5 \'^' ) 

4Lr5 

where ^trapped is the area of S. The condition that iS is a marginally trapped surface 
(meaning that it has zero expansion) can be translated into an unusual boundary problem: 
the functions need to satisfy 

3 



and 



V^, - ) (^± - $±) = (38a) 



= (38b) 

c 

= 4 (38c) 



where C is the curve (within the collision plane x~^ = x = 0) on which the surfaces 5+ and 
iS_ intersectn The indices a, b run over the direction, and g""^ is the inverse of the metric 



(34). 



Solving (38) exactly is difficult when the collision is off-center. However, an appropriate 



perturbative expansion leads to an analytically tractable problem. Consider the parameters 

C± = ' (39) 

and defin^I] 

By boosting in the x^ direction, we can change the ratio / without changing z+ or . 
In particular, we can choose a boost parameter such that E^ \Ej^ and E^ —>■ X^^E^ with 



A = E_/E^. After such a boost, we are in a frame where C+ = C- = C- We will generally 



^As before, we use the words "curve" and "surface" despite the fact that C has two dimensions and S has 
three. 

■^For comparison with [2 , it is useful to note that for central collisions in AdSs with z+ = z_ = L, 
= 4g(l + q){l + 2q) and x = 2yJ q{l + q). So when C is large, ( « 2g « a;. 
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prefer to work in such a frame. 

In p] it was noticed that in the case of central colhsions of identical objects, i.e., (3 = 
a = 0, the size of the trapped surface grows linearly in ( when ( is large. It is plausible that 
the same is true of non-central collisions, at least when the impact parameter f3 is held fixed 
while ( is made large. Let's define a new radial coordinate: 



P 



r 



(41) 



The metric on i/q takes the form 



dp' 



p2 + 1/(2 



+ CV(f/r + sin' 9d(j)' 



(42) 



We have already remarked that the off-center collision respects the 0(2) symmetry generated 
by additive shifts of (p. So the trapped surface must be a surface of revolution in the cj) 
direction. The curve C must likewise have the 0(2) symmetry, which means that its position 
in H-^ can be parameterized as p = pc{d)- Because we have chosen to work in a frame where 
Ej^z^ = E_z^, there is an additional Z2 reflection symmetry which interchanges the shocks 
by sending 6 ^ n — 6 and at the same time —x^. As consequences of this symmetry, 

we have 

^.{p, 9) = ^+{p, 7r-9) pc{9) = pc{n-9). (43) 



If we define 

Lh±ip,9) = ^±{p,9)-<!>±{p,9), 
then the functions h± are eigenfunctions of the laplacian on H^: 



L2 



p'dl + 3pd, - 3 



1 / 1 



C^p^ \sin^ 



-de sin 9de + p'^dp + 2p9^ 



The boundary conditions (38b) and (38c) can be re-expressed as 



h± + 



1 / 



L 

do ( h_ 



T 



(44) 



/i+ = 0, 



(45a) 



(45b) 



(45c) 
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The functions h± obey the same symmetry relation as h^{p, 0) = h^{p, n — 6). 
To make (45) analytically tractable, we expand 



h+{p,e) = hoip,9) + -h2{p,e) + ... 
Pcie) = poie) + ^p,{e) + .... 



(46) 



The differential equation (45a) can now be broken down order-by-order in (. At leading 
order in ( we find 

{p^dl + 3pd,-3)ho = 0. (47) 



The solution to (47) is 



hoip,9) = Coi9)p + Do{9)p~ 



(48) 



The second term on the right hand side of (48) can be discarded. Naively, the reason is 



that it is singular at p = 0. This is an unsatisfactory argument because (p = r/L, so that 
working in the large ( regime implies that r must be large. In order to impose the correct 



boundary conditions at small p we need, for example, to match the small p behavior of (48) 



to the large r asymptotics of the full solution to (45). We do this in appendix [A| where we 
find that the naive expectation Dq[9) = holds. Thus, 



ho = Co{9)p. 



(49) 



At subleading order in ( we have 



{p^dl + 3pdp - 3) h 



{p^dl + 2pdp + dl + cot 9de) ho . 



(50) 



The solution to (50) is 



h. 



C2{9)p + - [dl + cot 9de + 2) Co{9) . 
4p 



(51) 



As was the case for ho, the solution to (50) may also include a D2{9)p"^ term. In appendix 
|A]we also show that D2{9) = 0. One should be able, in principle, to carry out this procedure 
to arbitrary order in (. We find that at relative order some of the homogeneous solutions 
proportional to p^^ (which we may call D^ifi) following the notation in the preceding para- 
graphs) may not vanish. Additionally, there will also be a contribution at order C~^\ogC,. 



So going beyond the order shown explicitly in (51) seems to present some new difficulties 
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Fortunately, for practical purposes, ( is numerically fairly large: C ~ 50 or more in the cases 
we'll be considering. So we will not concern ourselves further with higher order corrections, 
and work only through relative order 

To determine the remaining integration constant Ci{6), we need to implement the bound- 
ary conditions (45b) and (45c). At leading order in ( |45b ) reads 



PoCo{9) + 



2pl (Vl+/32-/3cos^y 



0. 



(52) 



Plugging this into (45c) we get 



Po{0) 



1 



y/l + (3^ sin^ e 

This is the leading order expression for pc when ^ 1. At subleading order we find: 

1 -l + /?2(4cos(2^) -3) -6/5Sin2^ 



(53) 



Pc{e) 



+ sin^ e 



+ 



6(2(1 + /32sin2^)3/2 



+ 0{C 



(54) 



Reading off the volume element on from the metric (42), one can compute the area 
of the trapped surface S from 



A 



trapped 



2n 



" r^^'^ , l'Cp^ sine 

dO dp 
Jo 



(55) 



The explicit factor of 2 that multiplies the integral in (55) comes from the fact that S consists 



of two parts, iS+ and iS_, each with an area equal to the integral of the volume element 



over the region p < pc{0). Plugging (54) into (55) and using (37), we find that to leading 
order at large C, we have 



•S* ^ 'S'trappcd 



1/3 



TT- 



sinh ^ (3 



(56) 



An expression for (3 in terms of the impact parameter h = and the rms sizes z^ and 



Z- of the flattened energy distributions in the boundary theory (see (25) and (26)) can be 



obtained from (31): 



/3 



1 lh ^ + {z+-z_] 

Z+Z- 



(57) 



Plugging this expression into (56), one can obtain an expression for S'trappcd only in terms of 
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quantities defined in the boundary CFT. The special case of (56) where 2;+ = z_ 
= = E is equivalent to the result p|) that we quoted in the introduction. 



and 



Subleading corrections to (56) can be straightforwardly worked out from equations (54) 



and (55). One obtains 



trapped 



TC 



sinh (3 ^2 



logC - log2 



2 (1 + 2/32) sinh"^ /3 

3 ^ 2/VfT? 



{hi 



where C is as defined in (40). The non-analytic logC term in the above expression comes 



from the fact that acts as a small p regulator in the area integral (55). We can estimate 



that (58) will break down when /3 ~ because then the constant term in (58) becomes of 



the same order as the leading term, so one would expect that higher order terms in the series 
will then also be important. If, instead of colliding point-sourced shock waves in AdSs, we 
collide shock waves whose sources are spread out in the transverse plane, then we expect 
that the trapped surface will be almost unaffected provided it extends over a much bigger 
region of Hj, than the sources do. There is a precise result along these lines for head-on 
collisions [2]; see also the related discussion |10] . 

The area of the trapped surface was computed numerically in [36] for particular values 
of C, with = L. In figure [T] we compare the numerical results of p6] (in the case where 
G^E/L^ = 100, corresponding to C ~ 5.848) to the analytical prediction (58). As can be seen 



from this figure, equation (58) is a good approximation whenever 6/L < 4 (corresponding to 
/5 < 2), and breaks down for larger values of b. In particular, the approximate bound (58) 
doesn't capture the fact that there exists a maximum value b = fomax above which marginally 
trapped surfaces of the type considered above no longer exist. This doesn't necessarily imply 
that in collisions with impact parameters larger than b^ax there is no black hole formation, 
since there could be marginally trapped surfaces elsewhere. The existence of 6max is however 
suggestive of an upper limit on the impact parameter for black hole formation. 

The values C ^ 6 used in j3H] are smaller than the range we think can be compared 
most naturally to heavy-ion collisions in the range of energies attained at RHIC. As we shall 
explain in section |4| we prefer values of ( more than ten times bigger. For such large (, 



even the first correction to (56) is insignificant for the range of b/L we will be interested in. 



namely b/L no more than a few. 

As one approaches &max, the flat space analysis of shock-wave collisions described in 
[H] suggests that the slope of S'trapped(&) becomes infinite. This type of behavior implies 
that at 6 = 6max subleading terms in the series expansion (58) are important, and by the 
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Figure 1: (Color online.) Comparisons between the numerics of [36] and the analytic for- 
mula (58). The black dashed curve represents the leading term in (58); the solid red curve 
corresponds to the first two terms in (58); the dotted blue curve represents the expression 
(58), which is correct up to a term of order 0{l/('^)] the green dots represent the numerical 
evaluations used in figure 3 of |36]; lastly, the vertical green line marks the place where, 
according to [36], the maximum impact parameter h^^y^/L occurs. We thank S. Lin and 
E. Shuryak for providing us with the results of their numerical evaluations. 
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previous discussion, the maximum impact parameter leading to black hole formation satisfies 
^max ^ C-^- Because oc E^^'^, cl natural expectation is &max ~ 

^1/3 fQj. large energies. This 

scaling seems consistent with the observation of [36] that for C < 6 one has 6max ~ E"" 
with a ~ 0.37. Because b^g^^/^L ^ 0.91 according to [36] for ( = 5.848, and because the 
difference between the scaling with energy found there and the one we expect is fairly small, 
it is plausible that the estimate femax ~ is within a factor of 2 of the correct answer for ( 
larger than a few. 



4 Entropy production in a strongly coupled conformal 
field theory 

The motivation for our analysis is to gain some insight into entropy production in heavy- 
ion collisions. This implies that we need to somehow translate our AdSs trapped surface 
results into expectations for QCD. Any such translation is perilous, because what gravity 
calculations in AdSs are really dual to is a strongly coupled conformal field theory, and QCD 
is not such a theory. The perturbative behavior of QCD at high energies is surely relevant 
to the earliest stages of a heavy-ion collision, and confinement is obviously relevant for late 
stages. Our assumption in attempting to compare AdSs calculations with heavy ions is that 
there is an intermediate regime where QCD is fairly strongly coupled and fairly close to 
conformal, and that the dynamics of this intermediate regime is crucial to the production 
of entropy. Having stressed that these are assumptions which may fail to some degree, we 



will attempt in this section a comparison of the trapped surface calculation (56) to data. 
Later on in section |5] we attempt to make a quantitative estimate of how a failure of these 
assumptions will affect our results. 

In order to make such a comparison, we must first fix all the parameters appearing on 



the right hand sides of (56)-(57). We use 

LVG5 = 1.9, (59) 



which we obtain by comparing the equation of state of the SYM theory to that of lattice 
QCD, as in [2j. We must also relate E±, z±, b, and S to experimentally measured quantities. 
For E± we use the total beam energy, which in the case of RHIC collisions reads 

E^ = E^ = 197 X . (60) 
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The parameters z+ and z_ give the rms-energy-density-averaged transverse radius of the 
nuclei. A typical energy distribution for the nucleus is given by a Woods-Saxon profile which 
has an exponential fall-off. We've fixed z± to fit the rms-energy-density-averaged transverse 
radius of a gold nucleus resulting from a Woods-Saxon profile, 



4.3 fm . 



(61) 



See [2] for details. Combining (57) and (61) we obtain 

h 



(3 ^ 0.12 



fm 



(62) 



for a gold-gold collision. For ^snn = 200 GeV, combining (39), (40), (59), (60), and (61) 
gives ( ~ 77. As remarked in section [3| this is large enough that just the leading order 



estimate (56) can be used for b/L not too large. From here on, we will consider only this 
leading order approximation. 



Putting everything together, we can rewrite the entropy bound (56) as 



S > 35000 



\/snn 
200 GeV 



2/3 



sinh ^ (3 



(63) 



for gold-gold collisions. As explained in the introduction, the amount of entropy produced 
after a collision of two heavy ions may be inferred from a measurement of the total number 
of charged particles reaching the detector: 



S ^ 7.5N, 



ch - 



(64) 



Combining (63) and (64), we get 



A^chargcd > 4700 



2/3 



sinh ^ P 



200 GeV; 



(65) 



Most of the rest of this section is devoted to confronting the simple formula (65) with data. 



The impact parameter b can be related to the total number of nucleons participating 
in the collision, iVpart [12]. Our description here of how this is done in the Glauber model 
closely follows |13| HU |15]. Each heavy ion is replaced by a distribution of nucleons, which 
is proportional to the energy density and is given by a Woods-Saxon profile: see figure [2] 
We work with a distribution p proportional to energy density but normalized to unity. The 
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Figure 2: The energy density of a gold nucleus according to the Woods-Saxon profile, e = 

zr, n^sTT^: as a function of the radial distance r to the center of the nucleus. The energy 

distribution was normalized so that its value at r = is one. For a gold nucleus, the 

0.159 GeV/fm^ R = 6.38 fm, and a = 0.535 fm 



parameters R and a take the values eo 



effective distribution in the plane orthogonal to the beam axis (also normalized to unity) is 
given by 

T{r) = / p{z,r)dz, (66) 



where z is the direction of the beam and r* is a vector in the plane orthogonal to the beam 
axis. To determine A^part, one asks how many nucleons would scatter if there were no other 
interactions among nucleons other than inelastic nucleon-nucleon scattering, whose cross- 
section (Ttvat is measured. In the optical approximation, where one ignores the discreteness 



of individual nucleons in the distribution (66), the result is 



iVpart(&) = Na / T^(f) 1 - 1 - TB(f - b)ar,N 



+ Nb j fsif^ (l - - ^^(^"'- b)(rNN)j d\ , (67) 

where Na and A''^ are the numbers of nucleons in each nucleus. We took the values of cttvat 
from 



A detailed derivation of (67) can be found in |45j . 
As we have remarked, one usually uses a Woods-Saxon profile to obtain the transverse 



distribution (66) employed in Glauber model calculations. But the transverse energy distri- 



bution dual to a point-sourced shock wave is instead given by (25), which we will describe as 



conformal because it preserves an 0(3) subgroup of conformal transformations. It would be 
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b(fm) 




Figure 3: (Color online.) The impact parameter 6 as a function of the number of participat- 
ing nucleons A^part in a gold-gold collision, as obtained through optical Glauber calculations 
at a/s^v7v = 200 GeV, where a^N = 42 mb. The blue curve is based on the standard Woods- 
Saxon distribution, whereas the red curve is based on the conformal distribution, propor- 
tional to (25). Note that in going from 130 GeV to 200 GeV the scattering cross section 
decreases to 41 mb. 



more faithful to the AdSs computation to use the conformal distribution instead of Woods- 
Saxon. Perhaps surprisingly, for the energy ranges we are considering both these profiles 
give rather similar results. See figure [3] We therefore employed a Woods-Saxon profile for 
the Glauber calculations used to compare (65) to data. 

In figure |4] we compare our lower bound on the entropy as a function of the impact 
parameter with the PHOBOS data at ^/sNN = 200 GeV and ^/snn = 130 GeV, taken from 
[Hj (see also jlBJ). Instead of plotting Ai'charged in terms of the impact parameter b, it is more 
common to plot Ai'charged /Apart vcrsus Apart- Wc have done so in figure |5} 

The mismatch in b dependence between string theory results and the data is due to 
our not handling infrared effects correctly. The problem is that in a conformal theory, 
widely separated objects interact more strongly than they do in QCD. When two heavy-ions 
collide in the real world, the nucleons which are not in the collision region (spectators) are 
not expected to interact. On the other hand, in the collisions we have discussed, energy 
distributions which are far apart during the time of the collision will produce entropy. We 



can mimic the effect of spectators by setting and E_ in (60) equal to the fraction of the 
energy of each nucleus that participates in the collision according to a Glauber analysis. In 
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Figure 4: (Color online.) Total number of charged particles Ai'charged as a function of impact 
parameter b. The data was taken from the PHOBOS experiment |17]. The red curve 
corresponds to the lower bound on the number of charged particles which is based on the 



gauge-gravity duality (56). 
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Figure 5: (Color online.) Plots exhibiting the linear dependence of the total number of 
charged particles Ai'charged on the number of participating nucleons A^part- The data was 
taken from |17]. The shaded red region shows the allowed values of Ai'charged/^part, based on 
(|56|. 
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Figure 6: (Color online.) Plots exhibiting the linear dependence of the total number of 
charged particles Ai'charged on the number of participating nucleons A^part- The data was 



taken from [47J. The shaded red region shows the allowed values of A^charged/A^part, based on 
rescaling the right hand side (65) by f^^fe^) as explained in the text. 



other words, we rescale the energy of a collision with impact parameter b by 



Afpart(fe) 

2x197 ' 



where 



Npartib) is computed via the Glauber model. As a result, the right hand side of (65) is 



multiplied by 



2x197 J 



. The results of such a reinterpretation of the collision energy are 



shown in figure [6j 

While the fit exhibited in figure [6] is in good agreement with experiment, the approach on 
which it is based is somewhat ad hoc. A more consistent approach to fitting a holographic 
prediction with the data would be to repeat our calculation in a holographic dual of a 
confining gauge theory. 



5 Slicing trapped surfaces 

In the previous section we saw that the absence of confinement in our holographic model 
leads to a spectator problem in off-center collisions: too much energy far from the center of 
mass becomes thermalized. We solved this problem by redefining the energy of the shocks as 
the energy participating in the collision. We would prefer a solution where total energy can 
be cleanly identified between QCD and the five- dimensional bulk, and where the important 
non-conformal aspects of QCD are captured holographically. As a first step in this direction, 
we discard the parts of the trapped surface which lie outside a "safe" region 

ZiK> Z> ZuY . (68) 
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Our rationale for excluding the region z < zuv is that this region corresponds to the ultra- 
violet physics of QCD, which is asymptotically free. Partons are nearly free at high energy 
scales, so they often pass by one another without producing significant energy through scat- 
tering. Of course, there are hard scattering events, and these will be very interesting in 
LHC heavy-ion collisions. We are essentially ignoring such events, or folding them into the 
determination of the parameter zuy. 

We exclude the region z > zi^ because AdSs gets cut off in the infrared by confinement. 
The part of the trapped surface that has z > zin is not relevant to QCD because it relates 
to processes with energy lower than the confinement scale. Entropy production in real QCD 
must happen to some extent at energy scales below confinement, but it is presumably a small 
effect compared to the entropy production in the deconfined phase provided the total energy 
of the collision is large enough. 

In order to choose sensible values for and zyy, we first recall that the AdSs-Schwarzschild 
solution takes the form 

ds^ = ^ l^-hdt' + dx' + , (69) 

where h = 1 — z'^/zjj and L is the radius of AdSs. The temperature is T = I/ttzh- Based 
on this last relation, we associate physics at a scale A with the region of AdS^ with 

For the calculations in this section, we choose the values of z\jy and zib^ that correspond 



through (jrOj) to Auv = 2 GeV and Air = 0.2 GeV. 

In summary: We start with the trapped surfaces obtained in section [3] and intersect them 
with the "safe" region (68) to obtain a reduced entropy that excludes contributions from the 
ultraviolet, where we do not trust supergravity, and the infrared, where the bulk geometry 
should be cut off by effects dual to confinement. 

In the case of head-on collisions of objects of equal sizes, finding the area of the trapped 



surface inside the safe region (68) is fairly straightforward due to the 0(2) symmetry of the 
trapped surface in the x^x^-plane. Collisions which are not head on and/or those between 
objects of different sizes, are more challenging to handle. In the rest of this section we set 
z+ = z_, corresponding to equal-sized objects, although it should be possible to generalize 
our computations to unequal z^ and z_. We will work at an arbitrary impact parameter b. 
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We find it convenient to switch to a different coordinate system on defined through 



X-^ = Ljl+vl + vlcoshh 



X^ =L^l+vl+vlsmhh 
X^ = Lv, . 



(71) 



These coordinates have previously been used in 



When h, Vr, and v,, are unrestricted, 



they cover all of in a single coordinate patch. In fact, the mapping between and 



H'^ given in (71) is, up to a global rescaling, a volume-preserving diffeomorphism, as can be 
checked by noting that in the {h,Vx,Vy) coordinates the metric. 



1 + ^2+^2 



'1 + Vy)dvl - 2vxVydVxdVy + (1 + vDdv"! 



+ L\l + vl + vl)dh^ , (72) 



has determinant equal to L^. The 0(2) symmetry of the trapped surface manifests itself as 
rotations in the X2X^-plane, or as rotations in the fa;fy-plane. A constant h section of the 



trapped surface is a disk in the f^^fj^-plane. Keeping only the leading term in (54) and using 



the coordinate transformations (33) and (71), we find that the radius of this disk, v{h), is 
given by 



v{hf 



- 1 



1 + 2/32 + cosh(2/i) 
The area of the trapped surface, without cutoffs, is given by 



(73) 



A 



trapped 



dh Tcv{h)'^dh 



(74) 



To implement the UV and IR cutoffs of (68) we first need to find out what a surface of 



constant z looks like in the {h,Vx,Vy) coordinate system. Using (28) and (71) one can show 
that this surface is given by 



_ ^,center)2 (Vy - Vy 



center \ 2 



+ 



62 



(75) 
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Figure 7: (Color online.) A plot of a trapped surface (left) in the {h,Vx,Vy) coordinates for 
a head-on collision. The bright green circles correspond to surfaces of constant h. The blue 
surface (right) corresponds to a surface of constant z, whose constant h slices are ellipses 
shown in dark blue. At h = 0, the corresponding ellipse degenerates into a parabola. 



with 



center center N 



(0, (csch/i)2) 



(csch h)" 



6^ = (coth/i)^ 



(csch h)' 



for non-vanishing h. Thus, a constant h slice of a constant z surface is an ellipse whose 
eccentricity is 

(76) 



1 



cosh h 



When h = the ellipse degenerates into a parabola symmetric about the f ^^-axis, as can be 
seen from the fact that e = 1 and a = oo in this casej^ See figure [t) 

In figure [s] we've plotted the trapped surface (green) together with the z = ^uv surface 
(blue) and the z = zi^ surface (green). Pictorially, we need to compute the volume of the 



The parameters characterizing the parabola can be obtained from taking appropriate limits of (76 1 
the focal length is given by lim/j_^oO.(l — e) = and the focal point is located at w^. = and Vy = 
lim,,_o«™"'"' + ae = -2^-. 
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Figure 8: (color online) Plots of the trapped surface together with the UV and IR regions. 
On the left we show a plot of the trapped surface (green) in the {h, Vx, Vy) coordinate system. 
Some lines of constant h are lightly colored for emphasis. The UV surface and the IR surface 
are colored blue and red, respectively. Lines of constant /i 7^ on these surfaces are ellipses. 
On the right we've plotted a constant h section of the trapped surface and the UV and IR 
cutoff surfaces. 

green region which lies between the blue and red surfaces. In practice, it is easiest to compute 
the area of a constant h section T^h of the trapped surface restricted to zuv < z < zir and 
then integrate over h: 

POO 

A= Vol(S^)rf/i (77) 



where 

Sft = {{v^, Vy) : {Vx, Vy, h) eC and z^jy < z < Zyr}. {7i 
In figure [8] we show a typical constant-/;, slice. 



It is straightforward to see that the ellipse in (75) intersects the circle given in (73) at 



the Vy coordinate 



L V2(l+/32 + C2)cosh/i 



y™* = V ^ — (79) 

^ z ^1 + 2/32 + cosh 2/i 



The Vx coordinate of the intersection can be computed from either (75) or (73), but the 



resulting formula will not be needed. Let's denote by A]^ the two-dimensional area common 



to the disk of radius (73) and the region z > zir (which in figure [8] corresponds to the 
intersection between the disk and the interior of the inner ellipse) and by A]^^ the area 
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common to the same disk and the region z > Zuy (which in figure |8] corresponds to the 
intersection between the disk and the interior of the outer elhpse). The total area of the 
trapped surface is then 



Alice = 2 / dhAr - dhAi^\ . (80) 

\J — oo J —oo / 

The quantities A]^ and A]l^ can be computed analytically for every value of h. If, for 
example, the inner (IR) ellipse in figure [s] is contained completely inside the disk, then A]^ 
should just equal the area of the ellipse; if the disk is contained inside the ellipse, then A]^ 
should equal the area of the disk; if the inner ellipse intersects the boundary of the disk (as 
drawn in figure |8]), one needs to add the part of the area of the ellipse above the horizontal 
line passing through the intersection point in figure [8] to the part of the area of the disk 
below this intersection point. The explicit expression for AsUce is not very illuminating, so 
we won't reproduce it here. 

What we are interested in is how S'sUcc = ^siicc/4G5 compares with the entropy Strapped 
computed from the entire trapped surface. For central collisions, the results are shown in 
figure [oj converted to Ai'charged using the formula The UV cutoff doesn't significantly 

decrease the predicted bound on Ai'charged at RHIC energies, but as one proceeds to large E, 
its effect is dramatic: whereas Strapped oc E'^^^, we find SsHcc oc E^^^ in the limit of large E. 
The IR cutoff decreases the predicted bound on Ai'charged, but not by very much. To obtain 
the curves plotted in figure 9, we used L^/G^ = 1.9, z^: = 4.4 fm, and -Ebeam = 208 ^"^^^ ^ as 
appropriate for Pb. 



we 



Now let us turn to the case of non-zero impact parameter. In figures 10 and 11 
plotted the lower bound on Ai'chargcd corresponding to S > SsHcc for gold-gold collisions 
at y/SNN = 200 GeV and ^/sNN = 130 GeV together with the data from the PHOBOS 
experiment [121 SB]- The values for L^/G^, z±, and E used to make these plots are the ones 



given in the previous section. Figures [TO] and [TT] are an improvement relative to the fit in 
figures |4] and [5] in the sense that more data points are in the allowed region. Also, the string 
theory predictions at larger b are lowered by an incrementally larger factor than for small 
b. However, the string theory predictions still show too much entropy produced in glancing 



collisions as compared to central collisions. Comparing figure 10 to figure pi we conclude 



that our slicing procedure does not capture the effects of confinement as well as rescaling the 



^It should be noted that using (jlj over the large range of energy shown in figure [9] is not necessarily 
justified. Investigations along the lines of [55] based on LHC data should help clarify the relation between 
^charged and S at higher energies. 
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Figure 9: Total number of charged particles A'^charged for central Pb-Pb collisions over a range 
of energies. The red curve represents the 6 = limit of the AdS prediction (65). The 
blue curve represents the prediction of the Landau model (see, for example, section 2.3 of 
[2]). The brown dot-dashed curve represents the AdS prediction (80) with a UV cutoff at 
Auv = 2 GeV. The blue dashed curve represents the AdS prediction (80) with a UV cutoff 
at Auv = 2 GeV and an IR cutoff at Air = 0.2 GeV. The black vertical line marks the value 
of ^/snn expected to be attained at the LHC. 




Figure 10: Total number of charged particles Ai'charged as a function of impact parameter b. 
The data was taken from the PHOBOS experiment [121 HE]- The red curve corresponds to 
the lower bound on the number of charged particles which is based on the the dual of a slice 

of AdS dsel). 
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Figure 11: Plots exhibiting the hnear dependence of the total number of charged particles 
-^charged on the uumber of participating nucleons A^part- The data was taken from 
The shaded red region shows the allowed values of Ai'charged/^part computed via (80). 



energy did. It may be that a better treatment of the UV and IR cutoffs — in other words, 
better control over the non-conformal behavior of the theory — will yield a better justified fit 
to the data. 



6 Discussion 



Our main formal results, summarized in (54) and (58), are approximate analytic expressions 
for the size and shape of a trapped surface produced in off-center collisions of point-sourced 
shock waves in AdSs, in a limit where the trapped surface is much bigger than the radius of 
curvature of AdSs. The precise limit is C ^ 1 and ( ^ (3. It is quite striking that analytic 
results can be obtained, given that in flat space, the computation of trapped surfaces is a 
tricky numerical problem: see for example [¥T] . 

When using this formal result to make direct estimates of total multiplicities in heavy-ion 
collisions, we need to convert various AdSs quantities to equivalent QCD observables. If we 
interpret the energy of the colliding shocks as the total beam energy, then our results are in 
disagreement with the data at impact parameters greater than 4 to 5 fm. But if, instead, 
we interpret the energy of the colliding shocks as the energy participating in the collision as 
calculated via the Glauber model, then, as shown in flgure [6} agreement with data is good. 

Another challenge to our AdS-based model is the scaling S oc E'^^^ at high energies: the 
power 2/3 is probably too large. We've tried to address this problem by slicing off parts of 
AdS space. This crude approach is motivated by the idea that entropy production comes 
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predominantly from processes whose energy scale is above the confinement scale, but not 
so far above it as to be substantially suppressed by asymptotic freedom. Independent of 
the extent to which the slicing approach is justified, the calculations of section |5] provide 
information about where in the transverse H^, space (the slice of AdSs below the transverse 
x^x^-plane) the entropy comes from. 

To gain some intuition about the results of the slicing calculations, let's consider a head- 



on colhsion with 2;+ = z_ = L, = E_ = E, and ( = {2EG^/ L'^Y^^ , in accord with (39). 
First let's consider the shape of the trapped surface in the absence of any cutoffs. Let z^^-^ 
and Zmax be the minimum and maximum values of z on the trapped surface, and let x^^max be 
the maximum value of ■>/ {x^Y + (x^)^ on the trapped surface. Straightforward calculations 
to leading order in large C yield 

2^min 1 2^max > X_|_ j-q^x ^ \ 

~~r~ = — T — = ■ \p^) 



L 2C L 

If we compare to a central gold-gold collision at a/s/vw" = 200 GeV by setting L = 4.3 fm and 
L7G5 = 1.9, then ( = 77. (Note that setting L^/G^ = 1.9 makes C quite close to (EL)^/^ ) 



Combining the result (81) for z^-m with the conversion formula A = l/irz between depth and 
energy scale, one finds 

A„,ax = ^— ~2.2GeV. (82) 

In comparing to QCD processes, it is tempting to regard Amax as the maximum temperature 
of the fireball produced in a heavy-ion collision. But the trapped surface is far from being 
an equilibrated black hole horizon, so A^ax probably shouldn't be regarded as the peak tem- 
perature. Instead, it is the maximum energy scale of processes in the fireball that dominate 



entropy production. On this interpretation, the result (82) doesn't seem to us unreasonable 
in comparison with QCD. What seems more troublesome is the result for a;_L,max in (81): 
plugging in L = 4.3 fm and C = 77 gives 

a;±,max = 330 fm , (83) 

which, as noted in |36], seems way too high for a heavy- ion collision. 

Let's take a closer look at the distribution of the entropy in the transverse plane. For 
a central collision, it is fairly straightforward to calculate the fraction f{x±) of the entropy 
that lies within the region 

{x'f + {x'f < xl . (84) 
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It is interesting to note that although f{x±) is non-zero all the way out to x^^^^x = (L, its 
second moment is much smaller: 



/ /"Six, max 2 

x±,rms = Jj dx^xlf'ix^) ^ L^/2log2C - 3 ^ -logSEL - 3 , (85) 
where in the last expression we recalled that ( ^ [ELy/^ when L^/G^ = 1.9. We explain 



how to obtain (85) in Appendix pi With parameters as described above, x±^rms is roughly 



2.7L ^ 11 fm. Recalling that L is the energy-weighted rms radius of each shock wave, we 



see from (85) that the entropy-weighted rms radius x_L,rms is larger only by a modest factor, 
even at very large (. As we explain in Appendix |B| for x±/L of order unity, the distribution 
of entropy in the transverse plane as quantified by f{x±) is not so distant from sensible 



expectations for QCD. The surprising result (83), then, is not disastrous: it indicates that 



the tail of the entropy distribution is too long, not that the dominant part of the entropy 
distribution is unreasonable. 

Although f{x±) has a precise meaning on the gravity side of the duality, it is probably 
only approximately correct to translate it into the fraction of the entropy in the boundary 
theory within a radius x± of the collision point. Indeed, all results describing the position 
of the trapped surface in AdSs should be interpreted with caution when passing to the 
field theory side of the gauge-string duality. The trapped surfaces we have constructed 
stretch from t = — ootot = 0, but it would be obviously at odds with causality to say 
that the entropy is produced before the collision. In general, a black hole horizon is not 
something one can describe locally in time. It has to do with whether particles in a region 
of spacetime can escape to some asymptotic infinity. Trapped surfaces also have a non-local 



character: although the differential equation (45a) is local, the matching conditions (45b)- 



(45c) are not. The trapped surfaces we have constructed identify a region of spacetime from 
which a test particle cannot emerge once it enters — modulo some conjectural points relating 
to Cosmic Censorship. In particular, the test particle can be separated in the transverse 
(x^,x^) coordinates by as much as x_L^max and still fall inside the trapped surface. If we use 



(83), what this says is that a test particle more than 300 fm away from the colliding nuclei 
can eventually thermalize with the resulting fireball. Thermalization of such a distant test 
particle is impossible in QCD because of the short-range nature of the interactions. This 
underscores once again the real differences between real-world QCD and a strongly coupled 
conformal field theory. 

In Appendix [Bl we show that the IR cutoff lowers x_L,max from 330 fm to 13.8 fm in a 



31 



central collision with the choice of parameters indicated after (81). This dramatic reduction 
of a;_L,max is associated with less than a 20% reduction in the total entropy. In addition, 
a^±,rms is reduced by the cutoffs from 11 fm to roughly 5 fm for the parameters described just 



after (jSlf. 

In summary: The difficulty we had matching our results to total multiplicity of off-center 
collisions indicates that the trapped surface in the dual of a strongly coupled conformal field 
theory has too great a tendency to swallow up regions of spacetime that are far from the 
collision. In the conformal field theory, the fireball produced in an off-center collision has 
too great a tendency to thermalize with energy far from the collision point. Such energy 
corresponds to the spectators in an off-center heavy-ion collision. We suspect that this 
tendency goes hand in hand with excessively long tails of the entropy distribution in central 
collisions. Both difficulties would probably be cured in a holographic theory with a mass 
gap. 
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A Boundary conditions on h±{p, 6) at small p 



In this section we explain why we didn't include terms of the form Dij in (91). The 
starting point is the equation satisfied by which for simplicity we'll denote by h: 



Vl,3--^l/^ = 0. 



(86) 



Parameterizing i^s by (r, ^, 0) so that the metric is given by (34), the most general solution 



to (86) that is invariant under additive shifts in are 



h 



^>o 



f/(^)p^(cos 
+ \/(^)p^(cos0 



r 

T 



iFi 



- 1 



_+3 3 
2 '2' L2 



2F1 



1,1 



I 1 
2' 2 



r 

Z2 



^7) 



where P^(x) are the Legendre polynomials, 2-^1 is the hypergeometric function, and U^^'^ and 
are integration constants. Since the hypergeometric functions are regular at r = 0, it 
follows that the ?7- solution is regular at r = 0, while the V^-solution is not. The large r 
behavior is a bit more complicated: 



l>2 



U^'^PeicosO) (^a^^ + a?^ + 0(log r/r^)^ + vWp,(cos^) (^6^ ^ + 0{l/r' 



+ 
+ 

where 



v^r(2 - £) 
(£ + 2)r(i- 



^ 2^+iv^r(i - ^) 
,(.)_ (-irr(|-£)r(£ + 3) 



(89) 



As explained in section ([3]), equation (86) can also be solved in a coordinate system where 



Hj, is parameterized by [p^O^cf)) and the metric is (42). The parameter <^ = r/pL appearing 
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in (42) is a large expansion parameter. By expanding 



h{p,d)=h^{p,d) + -h2{p,d) + --- 



(90) 



one obtains 



10 



p^ 

1 1 



h2 = C2{e)p + D2{e)- + — (9,2 + cotede + 2) Co{e) 

p-^ Ap 



Up' 



(91a) 



(91b) 



The matched asymptotic expansion technique imphes that the small p limit of (91 ) needs 



to be matched onto the large r limit of (87) (equation (88)) after plugging in r = pL( in 



(88). To do so, we first write 



G 



2k 



cos( 



£>0 



D2k = J2D^^2Piicose) 



t>0 



(92) 



We find 



t/(o: 



a_i 

T 
1 

c 



^^2^2 + 
^2^2 + 



l(0) 



c 



2^2 



(0) 



hC [ D^'^ + -D'-> + ■ ■ ■ 



4(3 
3 
1 



D^o^ + i-2D\" + 



2CV ° ^e^^ ^ 



for £ > 2 



3 



^^2 2 



for £ > 2 



+ 4c^ ( /^r + ^^2" + 



(0) 



One can arrive at the first and fourth relations above by comparing the and 63"'' terms in 



(93) 



(88) to (91a) and the first two terms in (91b). The a\ term in (88) is the same as the third 



term in (91b) to leading order, but this is just a consistency check. A similar consistency 
check consists of comparing the fourth term in (91b) to the 0{l/r^) term that multiplies 



Actually, the next order term in (90 1 is of order C logC, so not all terms in this expansion take the 



form (911. 
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\/(^) which was not written exphcitly in (88). The other relations in (93) can be found by 



making similar comparisons of the terms in the second and third lines of (88) to (91) 



The main application of this matching procedure is that if one requires regularity of the 



0, one necessarily has Dq 



D 



for all i. To see this, note that 



solution at r 

regularity implies V^^"^ = 0. From the fourth and fifth relations in (93), this immediately 
implies that D^^^ = D^^ = for £ > 1. The same is true for i = since the D^^^ and D^^ 
terms in the last equation in (93) dominate over the Cg^'' terms at large (. Note that from 



the form of the formulas in (93) one may infer that D^^ = for all £ > 1 and all k, but it is 



not clear, for example, that D^^ should vanish. 

B The transverse size of the trapped surface 

In this appendix we consider the computation of the fraction f{x±) of the area of the trapped 
surface formed in a head-on collision that falls within the region (x^)^ + (x^)^ < x^. For 
simplicity we set = = L and = E_ = E. Throughout, we work to leading order in 
large C = {2EG^/L^f'\ 

The trapped surface has the familiar 0(3) symmetry that preserves the point of impact 
on the transverse space i^s. The region (x^)^ + (a;^)^ < x\ preserves an 0(2) subgroup of this 



symmetry. As in section [5| it is advantageous to use {h,Vx,Vy) coordinates in the presence 
of a residual 0(2) symmetry. But because the particular 0(2) we are now considering is a 
different subgroup of 0(3) from the one preserved in off-center calculations, we use a different 



identification from (71): 



= LJl + v?. + vlcoshh 



+ (xi)2 + (a;2^2 



1 + 



= Lv^ = L 



= Lv„ = L 



X 



X 



X* = -Ljl + vl + vlsmhh= - ( -1 + 



(94) 



The trapped surface covers the region of H3 with 



{x'Y + {x'^y + {xy < ei^ 



(95) 
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to leading order in (. Denoting v = y^v^ + Vy, we see that (95) may be recast as 



V < v{h) = W ggQ]^ }i _ tanh h . 



where h must be restricted to the range —hQ<h<hQ, with 



(96) 



ho = sinh ^ ( ^ \og2( . 



(97) 



Next, we observe that 



1 + f 2 + f 2 



L2 (X-l+X4)2 

This imphes that the inequahty (x^)^ + (x^)^ < is equivalent to 

V < — , or n < log — . 

" ^e^'' - xl/L^ - ^ L 



(9^ 



(99) 



Because the volume form on H^, is just Li^dh A dv^ A dvy, the fraction f{x^) can be computed 



as the ratio of the volume in the {h, v^, Vy) coordinate space satisfying both (96) and (99) to 



the volume satisfying only (96). At large (, one finds 



2L'^x± 



1 - 



(100) 



From (100) it follows that 



L2 



dx^xif\xA_) ^ 21og2C -3, 



(101) 



which gives the result quoted in (85). It is worth noting that f'{x±) falls off as l/x\ but is 



cut off at x_L,max = (L by the square root that appears in (100). 



To go further, let us define the distribution of entropy over the transverse plane as 

'j^'S'trapped 'S'trappcd n/ 



dx^dx^ 



2nx± 



(102) 
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Using (100), we find that at large (, 



I ^'^trapped 1 A ^ (103) 

Strapped ^^XW {L^+xlyV CL^- 



If we assume that 



J2 q 

>Jtrapped ^2 



where Qs is the saturation scale p] then we find 



Ql , (104) 



L2 + V 

which again holds only in the limit of large (. This formula shows that the saturation scale 
has approximately a l/i^L"^ + x^) dependence which is cut off at x± = x_L.max = (L, where 
Qs(a;±,max) = 0. It is striking that essentially the same spatial dependence, Qs oc 1/ {R'^+x\), 
has been advocated from the perspective of the color glass condensate [50]. A difference is 
that the parameter R of [50] is expected to grow logarithmically with energy, due to spread 
of the saturating region into the outer corona, whereas in our case, L is simply a constant. 



The slow growth of the rms size (85) of the trapped surface is due mostly to the lengthening 
of the power-law tail. In summary: transverse distributions reveal a qualitative, or even 
semi-quantitative, similarity between trapped surfaces in AdSs and the CGC approach; but 
the reasons for slow growth of transverse size with energy differ. 

In the CGC approach, the saturation scale at xj_ = is expected to grow as a slow 



power: oc where A ~ 0.28 (see for example [5T]). Using (105) and the scaling 
C ~ E^^'^, we find instead Ql oc E'^^'^ up to logarithmic corrections. This discrepancy 
probably owes to the fact that the QCD scaling Qg oc -E'^ is essentially a perturbative result, 
and hence distant from the strong coupling regime in which we work. Clearly, this is closely 
related to the expectation in CGC treatments that multiplicities will scale as a smaller 
power of E than the E'^^^ that we find. Thus we hold out some hope that an appropriate 
treatment of the ultraviolet physics, beyond the supergravity approximation, might lead to 
a more systematically accurate account of saturation physics as well as total multiphcity in 
a holographic framework. 

Including IR and UV cutoffs changes somewhat the story presented above. A surface of 



^^The saturation scale is a characteristic scale of transverse momenta of gluons in the early stages of a 
highly relativistic collision of nuclei. 
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Figure 12: The distribution of the trapped surface entropy in the transverse plane for 
^/sNN = 130 GeV, 200 GeV, and 500 GeV, in the presence of UV and IR cutoffs. We used 
L^/Gs = 1.9, L = 4.3 fm, Air = 0.2 GeV, and Auv = 2 GeV, as appropriate for a comparison 
to gold-gold collisions. 

constant z is described by 



1 , (106) 



and a careful account of the intersections of the IR and UV cutoffs with the trapped surface 



and the x± cutoff (99) yields 



L2 + Z^ 



in the limit of large (. Our preferred choice of zir, corresponding to Air = 0.2 GeV, reduces 



the value of x_L,max from 330 fm in absence of cutoffs to 13.8 fm with cutoffs. In figure [12 

^ '-'trappc 



d^ S 

we show the lower bound for the entropy density in the case of gold-gold collisions 
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at y^SNN = 130 GeV, 200 GeV, and 500 GeV. In making these plots we set L^/G^ = 1.9, 
L = 4.3 fm, Air = 0.2 GeV, and Auv = 2 GeV as discussed in sections |4] and [5] The IR 
and UV cutoffs aher somewhat the scahng ^^^'J^'^ ~ V(-^^ + ^1)^ ^^e values of ^/sNN 
considered above, but this scaling becomes more and more accurate in the intermediate x± 
region at larger (. At x± = 0, ^xriS'^ approaches a C-independent constant in the limit of 
large (. Assuming again that the entropy density is proportional to Ql, one can see that 
this scaling is significantly different from the result ~ C ~ -E^^^ in absence of cutoffs, and 
from the CGC scaling Qg ~ E^'^'^. What this shows is that changing the way we deal with 
UV physics does indeed affect significantly the scaling of Qg with energy at x± = 0. 

Our approach of entirely ignoring the ultraviolet region of AdSs (the region where z < 
z\jy) is probably too abrupt. Better would be to have some understanding of how entropy 
production weakens gradually as one departs from the regime of validity of supergravity. 
Such an understanding seems far out of reach starting from first principles in string theory. 
It may be that CGC results on the scaling of Qs provide some hints about the correct stringy 
dynamics of the ultraviolet region of the bulk in the true holographic dual of QCD. 
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